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Abstract 

The leading color-singlet contribution to the fragmentation function for a gluon to 
split into spin-triplet S-wave quarkonium is presented. In the case of charmonium, we 
find that this color-singlet term is always negligible compared to the leading color-octet 
contribution. 





The dominant mechanism for the production of heavy quarkonium at sufficiently large 
transverse momentum is parton fragmentation [|l|. In Ref. [^, we calculated the fragmenta¬ 
tion functions for gluons to split into S-wave quarkonium states to leading order in a* and 
to leading order in the nonrelativistic expansion. The fragmentation function for the 
state was given in analytic form, but only numerical results were presented for the ^Si state 
due to the lengthy and cumbersome expression for the fragmentation function. In this Brief 
Report, we present the detailed formula for the gluon fragmentation function into the ^Si 
state. 

A rigorous theory of the inclusive production of heavy quarkonium has been devel¬ 
oped recently by Bodwin, Braaten, and Lepage 0. This approach is based on the use of 
nonrelativistic QCD (NRQCD) to factor the production rate into short-distance coefficients 
that can be computed in perturbation theory and long-distance factors that are expressed 
as NRQCD matrix elements. Using this formalism, the fragmentation function for a gluon 

to split into a quarkonium state X with longitudinal momentum fraction z can be written 
as 

Dg^xiz,^) =Y,dn{z,^){0^) , ( 1 ) 

n 

where are local 4-fermion operators in NRQCD. The short-distance coefficients dn{z,fi) 
are independent of the quarkonium state X. For a fragmentation scale /i of order of the 
heavy quark mass uiq, they can be computed using perturbative theory in asi^mq). The 
dependence on the quarkonium state X appears in the long-distance factors (O^). The 
relative magnitude of the various matrix elements for a given state X can be estimated by 
how they scale with irtg and with the typical relative velocity v of the heavy quark inside 
the quarkonium. Thus, the factorization formula (|^ for the gluon fragmentation function is 
a double expansion in as and v. To determine the relative importance of the terms in this 
formula, one should take into account both the scaling in v of the matrix elements and the 
order in ag of their coefficients. 

We now consider the fragmentation function for a ^Si state which we denote by V. 
In the color-singlet model for quarkonium production |^, only a single term in the expansion 
(|^ for the fragmentation function is retained. In the notation of Ref. , the matrix element 
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in this term is ((9}^(^S'i)). This matrix element is proportional to the probability for the 
formation of the state V from a point-like QQ pair in a color-singlet ^Si state. The leading 
contribution to its short-distance coefficient d[ ^^\z) arises from the parton process g* —>■ 
QQgg and is of order af. Using the velocity-scaling rules of 0, the matrix element {OXi^Si)) 
is of order so the contribution to the gluon fragmentation function is of order The 

factorization approach reduces to the color-singlet model in the limit v ^ 0, since all other 
matrix elements in the expansion (|ID are higher order in v^. 

The charmonium and bottomonium systems are probably not sufficiently nonrelativis- 
tic that matrix elements that are suppressed by powers of can be completely neglected. 
In the case of a ^Si state, there is one matrix element that is suppressed by a single power of 
relative to but it also has a short-distance coefficient d{z) of order There 

are several “color-octet” matrix elements that are suppressed by two powers of including 
(C>^(^S'o)), ((9^(^S'i)), and {OXi^Pj))■ Of particular importance in gluon fragmentation is 
the matrix element {O^ (^5'i)), because it has a short-distance coefficient of order as due to 
the parton process g* —> QQ. All other matrix elements that are suppressed by have short- 
distance coefficients of order or higher. The matrix element {OXQSi)) is proportional 
to the probability for the formation of the state V (plus other particles) from a point-like 
QQ pair in a color-octet ^Si state. The corresponding contribution to the fragmentation 
function is of order compared to for the leading color-singlet term. The two fewer 
powers of can compensate for the suppression of the matrix element by Keeping only 
these two terms, the gluon fragmentation function for the ^Si state at the initial scale 2mQ 
can be written as 

D,^v{z,2mo) = d[’^'Hz,2mQ){OU'‘S,)) + 2mo){C>p=S0) ■ (2) 

To leading order in the coefficient of the color-singlet matrix element in can 
be deduced from the Feynman amplitude for g* —> QQgg, where the QQ pair is produced 
in a color-singlet state with vanishing relative momentum. The square of the amplitude 
for this process can be extracted from a calculation of the matrix element for e^e“ —rjjgg 
[^. There are several typographical errors that must be corrected in eq. (5) of Ref. In 
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the sixth term on the right side, the factor (1 — — x^) should be (1 — /i^ — X 2 )‘^. In the 

last term, (1 — /i — xiY should be {1 — jj? — xiY, (1 — /x — X 2 Y should be (1 — /i^ — 0 : 2 )^, 
{K- ■ -h_)^ should be (i^_ ■ -h_), and the overall sign of this last term should be changed from 
minus to plus. Having made these corrections, one can reproduce previous results for the 
energy distribution for 7 * —>■ tjjgg [^. The fragmentation function for the ^Si state can be 
calculated by the same method used for the ^S'o state in Ref. The calculation is rather 
involved and we present only the dual result: 


Xf"‘h.2mg) = £ dr 

X (fi{r,y) + gi{r,y) 


(l+r)/2 


dy 


{r+z^)/2z (1 - yY{y - r)2(|/2 - r)2 

1 + r — 2y y — r + \/y^ — t\ 

Tog- 


i=0 


2{y — r)^Jy‘^ — r y — r — y/y'^ — 


( 3 ) 


The integration variables are r = Attiq/s and y = p ■ q/s, where p and q are the 4-momenta 
of the quarkonium and the fragmenting gluon and s = q^. The functions fi and gi are 


/o(h?/) = r^(l-I-r)(3-I-12r-I-13r^) — 16r^(l-|-r)(l-|-3r)j/ 

- 2r(3 — 9r — 21r^-|-7r^)|/^ -|- 8r(4-|-3r-|-3r^)|/^ — 4r(9 — 3r — 4r^)x/^ 

- 16(1 + 3r + 3r^)y^ -f 8(6 + 7r)y^ - 32?/^ , (4) 

/i(H 2/) = —2r(l-|-5r-|-19r^-|-7r^)x/ -|- 96r‘^{l+r)y^ + 8(1 — 5r — 22r^ — 2r^)?/^ 

+ 16r(7 + 3r)j/^ - 8(5 + 7r)j/^ + 32?/^ (5) 

f 2 {'i^-,y) = r{l + 5r + 19r‘^ + 7r^) — A8r‘^{l + r)y — 4(1 — 5r — 22r^ — 2r^)x/^ 

- 8r(7 + 3r)x/^ + 4(5 + 7r)x/^ - 16y\ (6) 

5'o(H 2/) = 'r^(l — r)(3-|-24r-|-13r^) — 4r^(7 — 3r — 12r^)?/ — 2r^(17-|-22r — 7r^)x/^ 

-|- 4r^(13 -|- 5r — Qr^)y^ — 8r(l -|- 2r -|- 5r^ -|- 2r^)y^ — 8r(3 — Hr — 6r^)y^ 

+ 8(1 - 2r - 5r^)y^ , (7) 

5'i(h?/) = —2r^(l-I-r)(l — r)(l-I-7r)j/ -|- 8r^(l-|-3r)(l — 4r)x/^ 

-|- 4r(l -|- lOr -|- 57r^ -|- Ar^)y^ — 8r(l -|- 29r -|- 6r^)?/^ — 8(1 — 8r — 5r^)y^, (8) 

92 ir,y) = r^(l-|-r)(l — r)(l-|-7r) — 4r^(l-|-3r)(l — 4r)x/ 
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— 2r(l + lOr + 57r^ + Ar^)y‘^ + 4r(l + 29r + 6r‘^)y^ + 4(1 — 8r — 5r^)|/"^. (9) 

The integrals over r and ?/ in (|^) must be evaluated numerically to obtain the fragmentation 
function at the energy scale y = 2mQ. The fragmentation function at higher energy scales 
y is then obtained by Altarelli-Parisi evolution. 

The short-distance coefficient of the color-octet matrix element in (|^) can be cal¬ 
culated to leading order in as from the Feynman amplitude for g* QQ. The result is 

i 

„ X 7ra,(2mQ) 

4 {z,2mQ)= (5(1-^). (10) 

Q 

The radiative correction of order has also been calculated recently by Ma 0]. 

The relative importance of the color-singlet and color-octet contributions to the frag¬ 
mentation functions can be determined by integrating the initial fragmentation function (^ 
over z to get the fragmentation probability at the scale 2mQ : 

f\z Dg^v{z,2mQ) = (8.28 X 10-^)^ti^^(O}"(3^i)) 

JO mq 

+ (1.31 X . (11) 

'^Q 

The value of the color-singlet matrix element ((!l}^(^S'i)) can be determined from the elec¬ 
tronic width of the vector meson state. In the case of the charmonium states J/il) and 'ip', 
the matrix elements ((!l]^(^S'i)) are approximately 0.73 GeV^ and 0.11 GeV^, respectively. 
The most reliable determinations of the color-octet matrix elements come from recent data 
on prompt charmonium production from the GDF detector at the Tevatron. The color-octet 

contributions are necessary to explain the magnitude of the cross section for prompt J/'ip 
and ip' production at large transverse momentum, and they also explain the shape of the 
transverse momentum distribution [Q. The values of the matrix elements ((!l^(^S'i)) that 
are obtained by htting the observed cross sections are approximately 1.5 x 10“^ GeV^ for 
J/tp and 4.3 x 10“^ GeV^ for ip' p. Using the value as{2mc) = 0.26, we hnd that the two 
terms in the fragmentation probability ([Tl|) are approximately 3.2 x 10“® and 1.5 x 10““^ in 
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the case of J/ijj and 4.7 x 10“^ and 4.3 x 10“^ in the case of iIj'. The color-octet term is 
larger by about a factor of 50 for the J/ijj and larger by about a factor of 100 for the We 
conclude that the color-singlet term is always negligible compared to the color-octet term 
for charmonium. 

In the case of bottomonium, the color-singlet term in the gluon fragmentation function 
is probably also negligible. Because of the running of the coupling constant between the 
scales 2mc and 2mf,, the suppression of the color-singlet term by al decreases its relative 
importance by about a factor of 2. On the other hand, the suppression of the color-octet 
term by decreases its relative importance by about a factor of 10. The net effect is that 
the color-singlet term is still likely to be more than an order of magnitude smaller than the 
color-octet term. 

This work was supported in part by the U.S. Department of Energy, Division of High 
Energy Physics, under Grant DE-FG02-91-ER40684 and Grant DE-FG03-91ER40674. 
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